On contribution of three-body forces to Nd interaction 
at intermediate energies 
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Abstract 



' Available data on large-angle nucleon-deuteron elastic scattering Nd — > dN below 

the pion threshold give a signal for three-body forces. There is a problem of sepa- 
ration of possible subtle aspects of these forces from off-shell effects in two-nucleon 
potentials. By considering the main mechanisms of the process Nd — ► dN, we show 
qualitatively that in the quasi-binary reaction N + d —* (NN) + N with the final 
spin singlet NN-pair in the S-state the relative contribution of the 3N forces differs 
substantially from the elastic channel. It gives a new testing ground for the problem 
in question. 
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An existence of the three-nucleon forces (3NF) is not doubted both in the 
standard meson-exchange picture [1] and in the ciral perturbation theory [2]. 
Their strength and detail structure are still under discussion [3]. An often used 
model of a 3NF is the 2ir exchange in the form called the Tucson-Melbourne 
(TM-3NF) parametrization [1] (Fig.l). At present nucleon-nucleon (NN) forces 
underestimate binding energies of the light nuclei [4]. The 3NF allow one to fill 
in part the gaps. Another signal for the 3NF gives elastic Nd scattering below 
200 MeV. Recent experimental and theoretical investigations [5-9] show that 
reach set of spin observables in this process gives a real chance to study various 
aspects of 3NF effects. So, the measured cross section, which is underestimated 
at the scattering angles 9 cm = 60° — 180° by the data-equivalent modern 
NN forces, is excellently reproduced by three-body calculations with the TM- 
3NF. The deuteron analyzing power Ay is also well reproduced [10]. On the 
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other hand, the nucleoli analyzing power AP y is still under discussion [11], [12]. 
Furthermore, the same approach fails to explain tensor analyzing powers from 
the precise dp scattering data [10]. From this observation the authors of Ref. 
[10] conclude that there is deficiencies in the spin structure of the TM-3NF. 

However, one should note that off-shell properties of NN forces contribute 
simultaneously with the 3NF. The problem is that these ambiguous aspects of 
NN forces can be studied in the 3N or many nucleon systems only (see recent 
discussion in Ref. [12]). A noticeable NN input dependence was observed in 
Refs.[10,13]. In this connection the authors of Ref. [13] proposed to measure 
the zero point of the longitudinal asymmetry in the nd total cross section in 
order to get more clear signal for 3NF effects. This letter is another step in 
that direction. Namely, to separate more definitely the NN and 3NF effects 
in the elastic Nd scattering one should investigate supplementary processes at 
almost the same kinematics but with a different relative role of the NN and 
3NF. As shown here, a substantial different NN-to-3N ratio occurs in the pd- 
interaction with formation of the spin-singlet NN-pair in the final 1 S'o state. 
We propose here to study the reactions 

p + d — > (pp) + n (1) 

and n + d — > (nn) + p at large cm scattering angle of the secondary nucleon 
and low relative energy of two protons (pp) or neutrons (nn) E NN < 3 MeV, 
when the 1 S , state dominates. The reaction (1) was not yet investigated exper- 
imentally. Complete theoretical analysis of this reaction is possible at present 
at energies below the pion threshold (214 MeV) in the framework of rigorous 
3N-scattering approaches [14,15]. As shown by Faddeev calculation, at initial 
energies below about 200 MeV rescattering of higher order is very important. 
However, around 300 MeV the first two terms in the multiple scattering ex- 
pansion are sufficient to describe the total nd cross section [7] . Therefore, the 
first Born approximation can be used as a qualitative estimation near the pion 
threshold. Within this approximation one cannot find an exact contribution 
of the 3NF. Nevertheless, we can show qualitatively that the relative role of 
the NN and 3NF forces differs in this reaction considerably from that in the 
elastic Nd scattering. First, using isospin invariance, we show that the Born 
3NF amplitude (Fig. la) of the reaction with the singlet pn^So) pair is sup- 
pressed by the factor of | in respect to the 3NF amplitude with the deuteron. 
On the other hand, the Born term of the one nucleon exchange (ONE) mech- 
anism (Fig.lc), related to the NN- forces, is not affected by isospin factors. 
Second, the ONE contribution is modified considerably due to suppression of 
the higher orbital momenta I ^ in the final NN-system at low E^n whereas 
it is not so for the 3NF. In additional, an important modification of the single 
scattering (SS) mechanism (Fig. Id) occurs in the reaction (1) in comparison 
with the pd — > dp. 
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The dynamics of the reaction (1) is discussed here by analogy with known 
mechanisms of the backward elastic pd-scattering [16,17] (see Fig.l). As we 
show below, the ONE+SS sum dominates in the pd — * dp process at kinetic 
energy of the proton beam T p ~ 0.2—0.3 GeV. At higher energies T p = 0.4—1.0 
GeV the double pA^-scattering (Fig.le) with the A-excitation (A) gives the 
main contribution [16,18]. For the NN^ AN amplitude we use here the it + p 
exchange model [20] depicted in Fig. lb. This model describes the measured 
cross section of the pp— > pnn + reaction [19] in the A-region. Since the A 
mechanism is an important ingredient of the 3NF amplitude [21], we identify 
here the A contribution with the 3NF one. The A-isobar is considered as a 
stable baryon, that is appropriate below the pion production threshold. 

Let us discuss the isotopic spin factors for the Born 3NF term of the pd — > 
(np)s,tP amplitude depicted in Fig.la. The 2tt exchange mechanism contains 
two terms corresponding to different values of the total isospin of the inter- 
mediate meson- nucleon system, T = | and |: 



A(pd -> (pn) sP ) = A s T=l/2 + A s T=3/2 , (2) 
A(pd -> (pn) t p) = A l T=xj2 + A^ =3/2 , 

where the spin singlet (s) and triplet (t) states of the final pn-pair corre- 
spond to the isospin T pn = 1 and 0, respectively. The isospin structure of the 
amplitudes Afp is given by 

A^{pd^ipn)p) = (2T + 1)^/2(2^ + l)(T pn 0~|~)x 



x I 1 2 2 I I 2 2 T P™ I B s,t 



The Clebsh-Gordan coefficients and 6j-symbols are used here in standard no- 
tations. The dynamical factor Bt in Eq. (3) does not depend on z-projections 
of the isotopic spins. Assuming the spatial parts of the singlet and triplet wave 
functions of the pn pair to be the same, i.e. = B^, one can find from Eq.(3) 
the following ratios 



4 s 4 s 1 

= ^T=l/2 = /1 T=3/2 = _1 

At At o- V 1 ) 

^T=l/2 ^T=3/2 ° 

After substituting Eq. (4) into Eq. (2), one finds 

■ A(pd^ (pn) 3 p) 1 
Ul A(pd^(pn) t p) 3' l ° ' 
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We stress that owing to Eq.(4) the ratio (5) does not depend on the unknown 
relative phase nor on the ratio of the amplitudes Ai/ 2 and A 3 / 2 of the virtual 
process nN — > nN in Eq.(2). The result given by Eq.(5) is valid not only for 
the A-mechanism (Fig.le), as was found in Ref. [22]. In fact, all intermediate 
states of the meson- nucleon system both for the isotopic spin T=3/2 and 
T=l/2 are taken into account in Eq.(5) including the A and N* poles and 
the nN continuum. Obviously, the relation (5) is valid also for the sum of 
the diagrams in Fig. la with different combinations of the isovector mesons 
(7T7T, np, pp, . . . ), as well as for the reaction pd — > dN*. 

For the isoscalar meson exchange (u, 77, r]', ... ) we find the ratio Rf = 1. 
The same ratio is valid for the ONE mechanism, Rf NE — 1. It is impossible 
to write a definite isotopic factor for the SS-mechanism because in this case 
the s/t ratios are different for the isoscalar (r = 1) and isovector (r = |) NN- 
amplitudes, which are mixed with an unknown relative phase in the upper 
vertex of the diagram in Fig. Id. 

A detail formalism for the amplitude of the reaction pd — > (NN)N in the 
framework of the ONE+SS+A model can be derived from the pd — > dp for- 
malism of Refs. [16,18] . For this aim one should make the following substitution 
into the matrix elements 

|v? d >^ V^\^V >, (6) 

where \ipd > is the deuteron final state in the pd dp, m is the nucleon 
mass and | > is the scattering state of the final NN-system at the relative 
momentum k in the N + d — > (NN) + N reaction. Since the S-wave gives 
the main contribution to the singlet (iViV) s -state at E^n = — < 3 MeV, 
one should omit the D-component of the final deuteron state > in the 
pd — > dp formalism [16,18] when making the substitution (6). Thus, one has 
to insert into the upper vertex of the ONE diagram (Fig.le) the half-off-shell 
amplitude of iViV-scattering in the 1 S , state, t s (q,k). This amplitude, as a 
function of the off-shell momentum q, is very close in shape to the deuteron 
S-wave function in momentum space, u(q), and has a node at the point q ~ 0.4 
GeV/c. The node is caused by the short-range repulsion in the NN potential. 
A similar node available in the wave function u(q) can be connected to the 
null of the measured deuteron charge formfactor Gc(Q) at the transferred 
momentum Q ~ 4.5 fm _1 [23]. The node of u(q) was not yet observed directly 
in any reactions with the deuteron due to a large contribution of the deuteron 
D-state. An important feature of the reaction (1) is a possibility to display 
the node of the amplitude t s (q,k) straightforwardly in the cross section at 
T p = 600 - 700 MeV and 9 cm = 180° [24,25]. At initial energies 100-300 MeV 
this node makes the ONE contribution vanishing at 9 cm = 100° — 130°. 

The results of our calculations performed within the ONE+SS+A model with 
the Paris NN-potential are shown in Fig. 2 for the pd — * dp process and in Fig. 3 
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for the reaction (1). The model describes rather well the pd — > dp cross section 
at T p = 150 - 250 MeV and 6 cm > 120°. The A(=3NF) contribution strongly 
depends on the cutoff parameters A ntP in the NN^ AN amplitude. We use here 
the values A w = 0.6 GeV and A p = 0.7 GeV obtained from the fit the data on 
pp — > pnn + and pd — > dp [20,18]. The sum ONE+SS underestimates the cross 
section at 6 cm = 110° — 130°, but this discrepancy is eliminated by adding 
the A-contribution (Fig. 2a), as was observed in Refs.[6,21]. The calculated 
analyzing powers A v y , A y , and A yy are in only qualitative agreement with the 
data (Fig.2b-d). 

In contrast to the pd — > dp process, the influence of the A mechanism in the 
reaction (1) is rather weak in the cross section, but more pronounced in the an- 
alyzing powers (Fig. 3). At the minimum of the cross section, 9™ m = 120° — 140°, 
the role of 3NF increases due to i) vanishing of the ONE amplitude and ii) 
rather fast decreasing of the SS contribution (Fig. 3a). Outside of this region 
the A contribution to the reaction (1) is smaller than in the elastic pd scatter- 
ing owing to the isospin relations. Note, within the ONE+SS approximation 
the behaviour of the vector analyzing powers A v y and A y is considerably dif- 
ferent in the reaction (1) as compared to the pd — > dp process (Figs. 3b, d). 
The reason is a modified structure of the SS-mechanism. Indeed, only the 
pn-scattering at small angle contributes to the upper vertex of the SS mech- 
anism in the reaction (1) [24]. On the contrary, both the charge exchange 
process pn — > np and the pp elastic scattering at small angles contribute to 
the pd — > dp process [16]. The A mechanism taken into account in addition 
to the ONE+SS sum, changes the analyzing powers noticeably (Fig. 3c, d) 

In conclusion, the ONE+SS+A model allows one to understand qualitatively 
the main features of the pd — > dp observables at T p ~ 200 MeV. Within 
this model we found that the relative contribution of the 3NF in the reaction 
A" + d — > (NN)(}Sq) + A" differs considerably from the elastic pd scattering. 
A sizable modification of the analyzing powers is expected in the reaction 
(1) in comparison with the pd — > dp, in particular, due to the 3NF effects. 
Future experimental study of the reaction (1) near the pion threshold and 
rigorous three-body calculations, complementing recent analysis of the process 
pd — > dp, can gain more insight into the 3NF. 
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Fig. 1. Mechanisms of the reaction p + d — > (NN) + N: a) - the Born 3NF am- 
plitude of the pd — > (pn)p and pd — > processes, b) — ir + p exchange for the 
NN — ► iVA amplitude; (c) - one nucleon exchange, (d) - single scattering, (e) - 
A-isobar excitation 
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Fig. 2. The cm cross section (a) and analyzing powers Ay (b), Ay(d), Ay y (c) in 
the pd-elastic scattering in backward hemisphere at different initial energies of the 
proton T p and the deuteron Td (T p = \Td). The results of calculations with the 
different mechanisms are compared with the experimental data from Refs. [26](»), 
[27] (filled squares), and [10] (open triangles and circles): ONE (dashed line), SS 
(dotted), A (dashed-dotted), ONE+SS (full thick), ONE+SS+A (full thin) 
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Fig. 3. The same as in Fig.2, but for the reaction p + d — ► {pp) s + n at T p = 200 
MeV and relative energy of two protons E pp = 3 MeV versus the neutron scattering 
angle 0™ m 
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